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Introduction 


In the description of sources of electromagnetic radia- 
tion it is customary to regard them as located either at a 
finite or at an infinitely large distance ( "celestial" sources 
An actual source is regarded as celestial if the parallax of 
the source is close to the precision of its measurement. It 
is proposed that phasometry of the radiation permits the de- 
termination of the moments of passage of the radiation at var 
ious points of cosmic space at which are located stellite ob- 
servatories (SO) with recording instruments. This method is 
used to localize the source. The assumption of motion is ne- 
cessary to resolve the problems of ambiguity of phasometry, 
leading to ambiguity in the localization. Such a model is 
applicable, for example, to bursts of gamma-radiation, which 
occur as arae and short-lived increases in the intensity of 
gamma radiation[l] . The methods of localization of such 
sources is similar to the methods used in hyperbolic naviga- 
tion (the LORAN system) [2]. 

The author acknowledges Dr. of mathematical sciences I. 
V. Estulin and the participants of the Seminar of the Mathe- 
matical Section and the Section of Astrophysics and Radio- 
astronomy of the Institute of Space Research of the Academy 
of Science of the Soviet Union for their interest in his work 
and helpful comments. 



LOCALIZATION OF NON-STATIONARY SOURCES OF ELECTROMAGNETIC 


RADIATION WITH THE AID OF PHASOMETRY 
l. The Problem of Localization and Fundamental Relationships 
By localization is meant the determination of the three 
coordinates of a source or the determination of the direction 


to a celestial source. It may be assumed that the surface of 
the wave of radiation (that is, the geometric multiplicty of 
points in space at the same phase of radiation) appears either 
as a plane, normal to the direction to a celestial source, or 
as a sphere, at the center of which the source is located. 

For the localization we employ the measurements of the moments 
t^ of the passage of the wave at the points in space r^, given 
in a certain coordinate system OXYZ. The possibilities of lo- 
calization are defined by the collection n of the employed 
base points ?i(that is, by the number of SOs used in the ex- 
periment) . 


In the case of the recording of a plane wave the quanti- 
ties and t-t are related by the expressions: 


7\* i v + o(r i -O=o, i = 1, n- 1 . 


Cl 

V V 


( 1 , 1 ) 

( 1 . 2 ) 


where is the direction to the source in the coordinate sys- 
tem OXYZ; 

ty is the moment of passage of the wave at the origin 
of OXYZ: and 

c is the speed of transmission of the radiation (an as- 
terisk denotes the inner product) 

From (1.1)- (1.2) the measured t^ and the given r^ define 
the direction ly. The quantity t^ appears only formally and 
need not be determined. 

In the case of the recording of a shperical wave the 
quantities r^^ and t^ are related by the following expression: 

•Numbers in the margin indicate pagination in the foreign text 
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where is the radius-vector of the source in the system OXYZ; 
ty is the moment of the burst of radiation. 

Here the ty does not appear simply 1 formally and must be 

defined together with Fy~ // - 

The system (1.3) includes within itself the solutions 

of the systems (1.1) and (1.2) as asymptotes. However, if it 
is known beforehand that the source may be regarded as celes- 
tial, it is preferable to employ (1.1)-(1.2). Therefore for 
the construction of an algorithm for localization we will em- 
ploy both systems. 

Let us take the origin of the coordinate system OXYZ as 
at the point r Q , designate it as 0£ni, and employ the radius- 
vectors, defined by it, as p^, and pass to the time t, calcu- 
lated from the moment t Q : 

T f -t,-t # . L.J-fTi, (1.4) 

P Y = / V^o- ty’ty-t „• h* S > 

We introduce the matrix R, the i-th column of which is 
the vector p^, and the vector T, the i-th element of which is 
H, and we rewrite (1.1) in the form: 

/t*( v +CT = O , 

and (1.3) in the form: 

lp v l=-C' T v . 

|Pt “Py I * I I = C y 

Since (1.8) may be regarded as independent of (1.7), in 
following we will regard only (1.8) as defining p . 
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fflSSS* 


( 1 . 6 ) 

(1.7) 

(i.a) 

the 



In order to evaluate the precision of the localization 
we will take into account the error At of the measurement 
of t and the errors Ap-j^ of the given vectors p^. Since the 
Pi are defined in the coordinate system 0£nC, the origin of 
which is located at the point r 0> the error in the designation 
of r Q is excluded from the analysis. It is only necessary to 
regard this error as systematic in the case of a spherical 
wave. The importance of the errors will be evaluated in /6 

the linear approximations, with the following relationships, 
obtained by differentiation of (1.6) and (1.8), employed in 
the linerizat ion : 


p; A 


dr • 

j 



* 


(1.9) 

% 

(1 .10) 


o" aL. __ 6 ,. 7* 

Pi dp. > 

(1-7 )* d§r - c <S. . , 

\ C i C o) Q X . CU VJ ’ 


( 7.7 f Mx. = 6 v, i* 

\^i c o' dpj v 1 f 


where 


(i.n) 

(IU2) 


-4 


/i ~^Py Py’Pj I t i - i, U- 1 , 

4 = Py /p t . 


(I. I 2a) 

(uas) • 


and where 6-y is the Kroenecker delta function. 

However, besides (1.9)-(1.10) it is necessary to consider 
in linearization the differential relationship impled by (1.2): 


i; d i, =o. 


(I .13) 


It follows from (1.10) and (1.12) that the precision of 
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the localization of a celestial source In the linear approx- 
imation is only Influenced by the projection of the errors 
Ap ^ on 

T £ — ^ 

Ap ( V =£ v Ap if t = Tl-i, (!.!♦) 

and the precision of the localizat ic n of the \ource of a 
spherical wave by the projection of the errors on T i : 

^ — 94 

Ap. =£. a p, , i=r,~m . (i. is) 

As a measure of the quality of the statistical character 
of the errors of measurement we will employ the deviation 

of the errors Ax^ and the deviations o*y ^ nd 0 i\ Of the 

errors Ap^Y and Ap^ 1 , assuming these errors to be irremovable. 
We will assume the errors of measurement of the time of pas- 
sage of the wave Ax^ to be uncorrelated to the errors of the 
coordinates of the base points Ap d . The correlation among 
the Atj, and among the Apfy and the we wm issume to 

be the "worst" possible case, that 3s, maximizing the size 
of the elements ^(V) of the covariant matrices Vj ^ v $/j 5 

of the errors Aly/ X , AT Y/ . p of the localization of celestial 
sources and the covariant matrices V p/x> Vp/p the errors 

APy/ T , APy/p of the localization of the source of a spherical 
wave, denoting the size of the errors of measurement of the 
time and the given coordinates of the SOs, respectively. 


2. Case of n a 2 

The solution to (1.6.' in this case will be the simple 
vector ly » equal to: 


”P t /P, cy »« + it/KSinct , 


( 2 . 1 ) 



where 


( 2 . 2 ) 


COS Ct --Ct,/ Pl . 0 4CX4.T , 

and where k is an arbitrary vector orthogonal to p^: 

X* p t = 0. (2.3) 

In th celestial sphere the collection of possible solu- 
tions forms a circle of radius a with its center at the 
point P^/p* It follows from (1.9), (1.10), and (1.13), that 
only the projections and A£p/ T of the errors Aly/ T 

and Al Y/p may be defoined in the direction 3 , orthogonal 
to and lying in the plane passing through fy and 

s=V(* «P t )/IV(*«-p,)l. (2 * 4) 


The size of these projections will be equal to: 


-CAtj/pj Since, 

(2.5) 


/8 

A^r/p^^Ap^p.Sinoi, 

and their mean-square deviation: 

(2.6) 

<V c< VP. sino{ * 

(2.7) 

• 

= £T p' / Pi sinct • 

(2.8) 


The sum of the mean-square deviations will equal: 
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t 


(2.8a) 



a z /p < Si ” a 


where 



( 2 . 83 ) 


With the help of (2.8a) we may define the half-width of 
the ring in the celestial sphere (take, for example, triple 
Ofc) within the limits of which the source is located. Clear- 
ly, the error of localization grows without limit upon the 
approach of the direction to the source to the line of action 
of the vector p^. 

Let us now look at the localization of the source of a 
spherical wave. The multiplicity of solutions of ihe system 
(1.8) forms the surface of one of the branches of a hyper- 
boloid of revolution. One focus of the hyperboloid is loca- 
ted at the origin of the coordinate system 0£nC, and the 
other at the point p^. The major axis of the hyperboloid 
is equal to |ct|/2. For x>0 the surface of solutions is 
that branch of the hyperboloid which encloses the focus p 1 , 
and for x<0 the one which envelopes the focus located at the 
origin of the doordinates. If 1*0 the hyperboloid is trans- 
formed into a plane surface normal to p^ and passing through 
p^/2. The asymptotic surface of the hyperboloid is the cir- 
cular cone, formed by the direction of (2.1), with its vertex 
lying at the point p-j/2. In this way the asymptotic solution 
to (1.8) becomes the solution of (1.6). 

The consequences of the continuous solutions of (1.8) 
appear practically more suited to the solution of the problem /9 
which is the reverse of the problem of localization : the 

confirmation by means of system (1.8) of the coordinates of 
possible astronomical objects (for example, a portion of the 
region of the surface of the Sun) which might be the source 
of the radiation. From this follows the consideration, that 
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the source may be located not only on the surface of the hy- 
perboloid, but also in the neighborhood defined by the errors 
of measurement. For a definition of the possible statistical 
dispersion of the source from the surface of the hyperboloid, 
we turn our attention to the errors of measurement. According 
to (1.11), (1.12), the errors of localization Apy/ T and Apy/ p 
are related to Ax^ and Ap^l in the following expressions: 


^ P V/r " CAr i • 

(2.9) 

-K 1 i'y/p ^ Pi 

(2.10) 


These relationships only define the projections of the errors, 
a s 

p Y/t and Apy/p , on the perpendicular normal to the 
surface of the hyperboloid (1.8) at the point p^ : 



O/'V <.l- 


i ~ nr n i . 


(2. II) 


The mean-square deviations of these projections are equal to: 


- c ' ' sin k P, /x ) . 
Op' /2 ). 


(2.12) 

(2.13) 


where p-^ is the angle between anc ^ (that is, ^he P ar allactic 
displacement of the source arising from the transfer of the 
point of observation from the origin of the coordinate system 
to the point pi: 


p j =2arcsin(’l< 1 --<,|/ii), . n-i. 


(2.14) 


The sum of the mean-square deviations o p will be equal to: 


(T, - / 2 sin p 1/2 


(2.15) 
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where 


0_ =/ o " 4 + c*cr ) 

s ' p 4 / 

3. Case n=3 


i v</a 


(2.16) 


We will assume, that 

rank (R) = 2 (3.1) 

since for a rank of R equal to 1 the problem is reduced to 
that of part 2. We introduce the vector ic, equal to 

(3,2> 

with the help of which we rewrite (1.6) in the form 

7 . [ rf] (3.3) 




where [Ric] is a compound matrix, and [ c ^] is a compound 
vector, with the unknown element X. On the basis of (3*1) 
we may write: 

(3.4) 


— T * 


=-[**] 


cr I 

U 1 . 


Having inverted the matrix (see Appendix, p.l) we obtain: 

(3.5) 

_ 4- x / 1 / ~ 

1 ' 2 

We now find the value of X fromthe condition (1.2): 


*P 2 ]c r . + [o* h] CT )/x a . 


^ ^ x j/i - 


c^D^t 


where 




Thus , for 


c*f *D 1 t < i 


(3.6) 


(3.7) 


(3.8) 


the system (1.6) has two solutions symmetrically placed about 
the base plane passing through the points 0, pi,and 03 * 
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For 


( 3 . 9 ) 


/n 


c*f * D 




= l 


the system (1.6) has one solution, lying in the base plane. The 
ellipse (3*9) bounds the possible area of measurement of t. 

This ellipse is inscribed in the rectangle |ct 1 |^ p-^, I c t 2 1 ^ P2* 
Because of errors the measure quantity t may pass beyond 
the permissible region specific to the localized source near 
the base plane. In such a case it is necessary to place t 


within the boundary of the region corresponding to the sta- 
tistical fluctuations of the errors. 

It follows from (1.9), (1.10), and (1.13) that errors 
of localization are related to the errors of measurement in 
the following manner: 



1 




( 3 . 10 ; 


-M’y 


yp.~ ,7 * 




( 3 . 11 ) 


Obviously, errors of localization grow without limit as 
HytK+l . The quantity is equal to zero where k= 0, that 

is, when the vectors p^and P 2 are colinear, or when X^ is or- 
thogonal to k, In the first case the entire system of mea- 
surement degenerates, since the rank of R is equal to 1. In 
the second case X = 0 and the system has a single solution, 
lying in the base plane, but its precision in the direction 
orthogonal to the base plane is undefined. 

It is possible to show (see the Appendix p.2-p.6) that 
for the worst possible corrleation between At, and At ? and 

00 • * 

between Ap^Y and Ap n Y the mean-square deviation of the er- 
rors (3*10), (3»H) will be: 

°W = 7f * ma*|v(0 tf M r (3.12) 


9 


/12 


— 


a t/i = 1% max l fc v * fa 4 °pl PX 


(3.13) 


With the help of the triangular inequality it is possible to 
obtain the more satisfactory but larger values for and 

°i/ p : 


°l/r 6 c <'Vp ( 1 +( VpJ/ 51 " P sin ? ’ 


(3.14) 



4 C 



/P t +a ^/pJ/Sm ( 3sinv. 


(3.15) 


where 0 is the angle between and P 2 

4» is the angle between ly. and the base plane. 
Designating the smaller value between p^ and p 2 as P min > 
we obtain the following final value for the sum of the mean- 
square deviations of the errors of localization: 


°t 4a 2 /, P"m’ sin P sit ' H ' - 


(3.15a; 



i ) 1/a * *(«v+ 0p) ,/a 


(3,156) 


where a T and cr are the larger of o o - and o , and o ,,, 

T P xl t 2 pi p2* 

respectively . 

We move now to the localization of the source of a spher- 
ical wave. According to (1.8), the multiplicity of possible 
solutions forms a curve defined by the line of intersection 
of two hyperboloids. This curve has two asymptotes, colinear 
to (3*5). As with n=2, it is approrpiate here to solve the /13 
reverse problem. For this, in order to define the area with- 
in the limits of which the source may be found, we turn our 
attention to the errors of measurement. 

The expressions (1.12), (1.12) relate the errors of lo- 
calization in the perpendicular (2.11) and the curve (1.8) 
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at the point solely to the projections of the errors of 

measurement and A Py/p* the vectors (2.11) are not 

linearly independent, that is if the source is not located 

along one side or the extension of a side of the triangle 

formed by the base points 0,p 1 ,p 2 , then it is possible to 

define the projection Ap p of the errors of the localization 

of the source on the plane passing through the normal vectors 

(2.11). Assuming, as before, the worst possible correlation 

o , o 

between At-^ and At 2 and between Ap.^1 and Ap g, it is possible 
to obtain the following values for the mean-square devia- 
tion of this projection Ap p : 

<y p /T £ c (< 7/2 sin (A/ 2 )+<j /2 sin( p V 2 ))/sitt Y , <3.i6) 
ff PA 4 ( <y p/~ s ' n ( Pl /3.) + <5pl /a sir. ( p Va ))/sin y , 

where y is the angle between X^-Z 0 and X 2 -E Q and p i ls def i ned 
in correspondence with (2.14). 


Analogous to (3.15a) and (3. 15b) we obtain the following 
values for the sum of maen-square deviations of the errors of 
locaxization : 


of* 0 2 /2»unp m , n /a) sin V , 


(3.1b) 


and for small parallax: 




trim 


Sin y . 


(3.19) 


4. Case n=4 


We will assume that: 


S5S5® 


li 



rank (R) = 3 


(4.1) 


that is, the base points do not lie in one plane. For the 
localization of celestial sources It is possible to employ 
algorithm (3»5)-(3*7), selecting the three of the four 
points which provide the most precision of localization, 
while using the remaining point to resolve ambiguity. How- 
ever, it is more satisfactory to immediately obtain the sole 
solution to the system (1.6), using condition (4.1). We in- 
troduce the vector S, equal to: 

s =-R*' l cr . ( 4 . 2 ) 

In the absence of errors of measurement the vector 3 would be 
the solution to (1.6). The influence of errors makes appar- 
ent the consequences of the redundance of (1.6) for n=4 and 
introduces a violation of the condition (1.2). Therefore the 
solution to (1.6) will be the normalized vector S: 

f v =V-S. (4.3) 

Expression (4.3), essentially, averages the measurements, in- 
cluding the redundance of the fourth measurement. Therefore 
the precision of the localization in this case can be greater 
than with the use of ( 3 • 5)-( 3* 7) • However, this greater pre- 
cision does not qualitatively change the preceding analysis, 
and the values obtained there may be emplayed in this case. 

Let us look at a more general algorithm, which includes 
the localization of both a celestial source and the source 
of a spherical wave. We square both sides of equation (1.8) 
and rewrite it in the form: 

g-p Y cf • (*.4) 
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where the elements of the vector q are: 


« T^rr. 


(4.5) 


Since |ct^| < pi , then: 

P*/z » ? 4 > c . 


14.6) 


/ TTTT 

As a consequence of (4.1) we may solve (4.4) with respect 

P v =P+pJ ’ 


to p y : 


Where S Is as defined in (4.2), but : 


P - * “ 1 9 . 

If all q^ equal zer .•>, then: 


(4.3) 


(4.9) 


and since according to this S / 1, the only solution to equa- 
tion (4.9) would be p^ = 0. On the other hand, from p^ * ft 
it follows that q = (3. Excluding from further analysis the 
case p^* C let us assume that: 

l cr i l < Pi > p*/a>y i > f. p y >o. (4.io) 


and let us introduce the variable: 





with the help of which we rewrite (4.7) in the form: 


I 



(4.12) 


where I Q is the single vector of the direction to the source 
from the origin of the coordinate system 0£nC (1.12b). 





OF 


If 3 = 5, that is if 7 = 13, then it follows from (4.10) 
that p 0 is equal to the vector p 0 : 



(4.13) 


Since for ct = 0 the source is located at an equal distance 
from all base points, the p Q is the radius-vector of the cen- 
ter of the sphere which passes through the base points 3, pj, 
P2 > P 3 * 


In the general case we obtain from (4.12) the following 
equation: / l6 

/> Vy 4/5 * s n- v + s* - i =o, p*c. (4#tt) 


For S<1 (4.2) has the single positive solution: 

fi. =(-p , 'S + (p‘ 4 - |/>x5 1' 3 ) 7 )/p a . (4.l c : 


In conjunction with (4.2) the condition S<1 satisfies the 
collection of vectors ct lying within the region of the 
ellipsoid : 


c*i *d l i - 1 . 


(4.16) 


For S^l both solutions of (4.14) may be positive. In par- 
ticular, when S = 1 (that is, when ct lies on the surface 
(4.16)), one of the solutions corresponds to a celestial 
source, the direction to which is : 


and the second, to a source located at the distance: 


(4.17) 


14 



( 4 . 18 ) 



For S>1 the solution lies on the various sides of the three- 
sided figure with its vertex at the point 3 and edges p^, P2, 
p^. For the solution lying on the interior side (that is, 
on the side on which is located the center of the sphere de- 
scribed by (4.13)), the quantity y^ is equal to ( 4 . 15 ) and 
0 < p y < p Q . For rhe solution lying on the exterior side of 
the figure, the value of y^ is equal to 

H- v =(-p “s-(p*-\p *5| a )‘ / ’*)/ p* 

and 0 < p^ < 00 . 


We turn now to the precision of the localization of the / 17 
source of a spherical wave. In accordance with (1.11) and 
(1.12) the errors Ap and Ap Y//p will be equal to: 

* *-4 

Ap v/t .=Z CAT, (4.20) 


• A p y jp /• A p , 


($. 21 ) 


wher L is the matrix whose i-th column is 

-I ^ £4 

and Ap is the vector whose i-th element is Apj. 1 

Thus, the system of measurement degenerates for 


(1.15). 
I L | = 0. 


The determinant |L| is equal to zero when the ends of the 
vectors Z Q , 1^, Z^ ,1^, constructed from one axis, all lie 
in the same plane. In particular, |L| = 0 if the source is 
situated along one of the edges or the extension of an edge 
of the pyramid whose vertices are the base points. In the 
general case only the quantities of the parallaxes p i (2.14) 
influence the precision of the localization, since the ma- 
trix L, and consequently the characteristics of the precision 
of localization as well, remain unchanged under the arbi- 
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ORIGINAL 1‘ • '2 .S 
Ofi POOR QUALITY 

trary rearrangement of the base points from any original ar- 
rangement along a straight line passing through the original 
base points and the source. 


In the case where the distance to the source is compara- 
ble to the distances between the base points, there is an 
optimal distribution of base points which maximizes the pre- 
cision of the localization, where the columns of the matrix 

r 

L are orthogonal and the parallaxes p i (2.14) are equal to 
2 arcsin (1//3). In this case, assuming the worst possible 
correlation (p.2-p.6), the mean-square deviation of the er- 
rors of localization Apy/ X and APy/ p are bounded by the 
following quantities: 


(T 4 e <T«. 
P/r a t * 


P/P 


C— G 
* a u p 


( 4 . 22 ) 

( 4 . 23 ) 


In the case where the distance to the source is significantly / 18 
greater than the distances between base points, the precision 
of localization is reduced proportionate to the square of the 
smallest parallax of the p 1 (2.14) (let it be p 1 ): 

***<**//>?. 



( 4 . 25 ) 


The sum of the mean-square deviations of the errors of lo- 
calization will be bounded by the following quantities: 


where 



( 4 . 26 ) 


_ , a 2 1 1 

V rc< V <J p, 



( 4 . 27 ) 
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Since p^,- p^/ p^, then for the relative error of localization 
we obtain the following expression: 

s C p/Py Pi (~P*t) • (k.2b) 

Assuming that the relative error may not, in any case exceed 
unity, we may conclude that: 

Pi ~Pi 4 K ’ K>i • (4.29) 


Thus, for a given precision there is a limiting distance 
p ymax* f ' or which we may ot) tain the statistically probable 
value : 


P Y way / 2 A # 


(4.30) 


5. Conclusion 

For the localization of a celestial source it is suffi- 
cient to register the burst at three SOs . Ambiguities which 
arise in this case may be removed by the use of directional 
detectors. With the use of detectors with an angle of di- 
rectional reception equal to 211 steradians, for unambiguous 
localization there must be in one SO a plane surface, co- 
inciding with the base plane, dividing it into two direc- 
tional detectors. Another possible means to resolve ambi- 
guity is to take steps leading to a fourth SO. 

Using three SOs it is impossible to distinguish a plane 
wave from a spherical one (this is only possible if the source 
is located on the base plane). Thus the treatment of the de- 
tection of a pulse from a nearby source defines a direction 
to a celestial source which may be altogether different form 
the direction to the true source. Therefore, in the use 
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of three SOs it is necessary to verify (the reverse problem) 
the possibility of radiation pulses from "suspicious" areas. 

If these regions are sufficiently small in one measurement 
(for example, a thin layer of the surface of the Sun), then 
in the process of solving the reverse problem one may define 
all three coordinates of the source. 

Id ntification of the form of the wave and the resulting 
definition of the direction to a celestial source or the 
three coordinates of a nearby source is possible with the use 
of four SOs. 

In the general case with the use of four SOs the coordi- 
nates of a nearby source are defined ambiguously. 

In all cases localization suffers a degeneration of the 
system of measurement if the source is located on a line 
extending through any two of the base points. In the treat- 
ment one of these points must be discarded, and localization /20 
based on the remaining points. 
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APPENDIX 


We introduce several relationships used in the funda- 
mental parts of the work. 

a. ) For the nondegenerate matrix A, its inverse may be 
represented as: 

A i 9 C ) * [I xc cxa 5 x £ ] ^ (p.l) 

where a, 5, c are columns in the matrix A*. 

b. ) Consider matrices of the type V = AKA*, where A is 
the given matrix and K is an arbitrary correlation matrix: 


K • = 

1 n 


*j 


r A 


jf 


I 




The largest result S of the matrix V is reached when K = K , 

max* 

where 


*Wij =s 'gn( a i Qj), ij’i.n, p- 2 > 

and equal to 


W S “P s(v) |q‘ a. 1 4 (f)|a t |') , 

A O t 


( P - 3 ) 


where a i 
The 


is the i-th column of A. 

maximum of the values of the matrix V is bounded by 


max 


p,<aa:a*)<. 

A 




max • 


(p. *0 


In particular, 


if the rank of V 


max 


u ^ S 

i ma-jf. 


max max 


is equal to unity. 


then : 

(p.5) 
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For n a 2 the rank of K max is always equal to one. For n>2 

u and S are constructed of various K „ . For example, 
M max max max * 

if n=3 there exist four different K „ . dependent on A, from 

max 

which S may be constructed: 
max 


'1 1 i “ 

i 

1 

i — 

i 


' 1 -4 -i' 

1 4 4 

» Aj - 

- 1 i ! 

> K„= ri i i 

. K 4 = 

“4 i - , 

i 1 4 

! 

1 

i 

8 L-* i *J 


.-1-4 i. 


(p.6) 


The rank of and K-, is equal to 1, but the rank of K.~, and 
is equal to 3* 
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